Disorder instability of the magnon condensate in a frustrated spin ladder. 
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The effect of disorder is studied on the field-induced quantum phase transition in the frustrated 
spin-ladder compound H8C4S02-Cu2(Cli-ccBrcc)4 using bulk magnetic and thermodynamic mea- 
surements. The parent material (x = 0) is a quantum spin liquid, which in applied fields is known 
to form a magnon condensate with long-range helimagnetic order. We show that bond random- 
ness introduced by a chemical substitution on the non-magnetic halogene site destroys this phase 
transition at very low concentrations, already for x = 0.01. The extreme fragility of the magnon 
condensate is attributed to random frustration in the incommensurate state. 

PACS numbers: 75.10. Jm,75.40.Cx,75.50.Lk,75.10.Pq 



I. INTRODUCTION 

Three-dimensional (3D) long-range order is robust and 
usually highly resistant to weak spatially random per- 
turbations. In particular, this applies to field-induced 
magnetic order in gapped quantum spin systems. This 
type of quantum phase transition has recently attracted 
a great deal of attention due to its interpretation in terms 
of a Bose-Einstein condensation (BEC) of magnons.^ ^ In 
magnetic materials, disorder can be introduced by ran- 
domizing the strength of magnetic bonds. In simple mod- 
els, its primary effect is to produce a random potential 
for the condensing magnons. The result is a qualitatively 
new phase with non-zero magnetization and susceptibil- 
ity, but only short range correlations, the so-called mag- 
netic Bose glass. ^ ^ However, unless the disorder is very 
strong, the ordered BEC phase is expected to re-emerge 
at higher fields.^ An intriguing question is whether weak 
randomness in quantum magnets can go beyond creating 
a random magnon potential, in a way that disrupts the 
formation of the condensate altogether? 

In the present work we investigate the effect of 
bond strength disorder on the geometrically frustrated 
gapped quantum magnet Il8C4S02-Cu2Cl4 (SCC for 
short) !^ Bond disorder is introduced by a chemical 
substitution on the non-magnetic site in the isostruc- 
tural H8C4S02-Cu2(Cli_a;Bra;)4 (hereafter abbreviated 
as sex). We find that extremely low Br concentra- 
tions have a catastrophic effect on the phase transition 
and magnetic properties. This behavior is attributed to 
frustration disorder^ enhanced by the incommensurate 
and quasi-one-dimensional spin correlations in the par- 
ent compound. 

To our knowledge, the parent compound SCC is quite 
unique being a strongly quasi-one-dimensional gapped 
Heisenberg spin system where geometric frustration is 
strong enough to cause dynamic spin correlations to be 
incommensurate. From the point of view of structural 
chemistry, it belongs to a large and diverse family of or- 
ganic complexes involving transition metal halogenides 
(see, for example, Refs. i9r rll ) . Antiferromagnetic (AF) 



interactions between the S = 1/2 Cu^+ ions are car- 
ried by Cu-X-X-Cu superexchange pathways (X=halogen 
ion). The resulting spin network can be described as 
four-leg "spin tubes" with a high degree of geomet- 
ric frustration,!^^ and weak inter-tube coupling. The 
ground state is a spin liquid, with a gap A ^ 0.5 meV 
in the magnetic excitation spectrum. The exact spin 
Hamiltonian is presently unknown. However, due to the 
strong one-dimensionality, all relevant low-energy low- 
temperature properties are governed by the lowest-energy 
triplet of magnons. in SCC, the latter have a very 
large spin velocity v along the c axis, hvc ~ 14 meV. 
A key observation is that, due to the geometric frustra- 
tion, the minimum of dispersion is at an incommensu- 
rate wave vector qoc/{27T) = 0.48. In magnetic fields 
exceeding He = 3.75 T a condensation of magnon^ 
produces an incommensurate phase with spontaneous 
3D long-range helimagnetic order of spin components 
transverse to th e app lied fields and a propagation vec- 
tor (0.78, 0, 0.48).™3This pitch of the spin spiral corre- 
sponds to the incommensurate dynamic spin correlations 
in zero field. 



II. EXPERIMENTAL PROCEDURES AND 
RESULTS 

To investigate the effect of disorder on this field- 
induced phase transition, we focused on the chemically 
modified material SCX, with the Br concentration x 
ranging 0-5%. The idea is that substituting Cl~ by the 
larger-radius Br~ will affect the Cu-Cu superexchange 
coupling locally.'^ Since the relevant Cu-X-X-Cu bond 
angles are far from 90°, Kanamori-Goodenough ruleJ^^ 
suggest that the AF nature of the interaction will be 
preserved, although its magnitude may be significantly 
affected. A random distribution of Br substitution sites 
will thus result in a randomization of AF bond strengths. 

A series of well-faceted single crystal samples of typ- 
ical mass 100 mg were grown by temperature gradi- 
ent method from ethanol solution with varying relative 
CI/Br content. For each nominal Br concentration x. 
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FIG. 1: (Color online) Connected symbols: specific heat measured as a function of temperature at various applied fields in 
sex with X = (a), X = 0.01 (b) and x — 0.035 (c). The heavy solid curves are guides to the eye emphasizing curves peaked 
at T 1 K for easy visual comparison. 



the samples were characterized by X-ray diffraction us- 
ing a Bruker Apex-II single crystal diffractometer and by 
micro-elemental chemical analysis (Schoniger method). 
We found all crystals to be isostructural to the parent 
compound. Atomic positions were refined together with 
the actual Br concentration, assuming a uniform distri- 
bution of substitution sites. The fitted Br content is in 
good agreement with chemical data and the nominal Br 
concentration in solution. As far as the lattice parame- 
ters are concerned, the strongest effect of Br substitution 
is on the c (ladder-leg) axis, which shows a linear increase 
from 10.035 A to 10.050 A, as x varies from to 5%. 



Characterization of the disorder-free material 

The clearest signature of the field-induced ordering 
transition in disorder-free SCO is the lamb da- anomaly 
in the temperature dependence of specific heat.^^ Fig.jTp 
shows the heat capacity measured in a 3.5 mg x = 
sample as a function of temperature in several fields ap- 
plied along the crystallographic c axis, using a Quantum 
Design PPMS with a dilution refrigerator insert. The 
sharp lambda anomaly is quite prominent and allows to 
trace the phase boundary on the H — T phase diagram, 
as shown in solid symbols in Fig. [2] The phase tran- 
sition is also manifest in the field dependence of mag- 
netic susceptibility x = dM/dH (Fig. [sj solid symbols), 
measured for H\\c on a 25 mg crystal using a Quan- 
tum Design SQUID magnetometer with an IQuantum 



^He refrigeratoi'^. The transition appears as a sharp 
peak, separating the gapped phase at low fields and 
the ordered, susceptible high-field state. For a BEG of 
magnons, magnetization M is mapped onto Boson den- 
sity, and susceptibility x - the compressibility of the 
Bose gas.^ In these terms, the transition is that between 
an incompressible "Mott-insulator" magnon phase be- 
low He, to a compressible Bose condensate of magnons 
at high fields. Note that an elevated temperature of 
800 mK, the step in the x{H) curve is broadened, but 
a sharp peak marking the phase transition persists. The 
third manifestation of the 3D ordering transition is found 
in M(T) curves measured at fixed fields applied along the 
c axis (Fig. |3j inset). Upon cooling through the ordering 
temperature, the magnetization (Boson density) shows 
an almost discontinuous jump following a well-defined 
minimum, which is a known characteristic signature of 
BEC.^ 



Effect of disorder 

The central result of this work is that Br substitu- 
tion at very low concentrations destroys the thermody- 
namic field-induced transition described above. Already 
for X = 1% the lambda anomalies in the C{T) curves 
at all fields are replaced by visibly broadened peaks 
(Fig. |T]3). For X = 3.5% the peak widths increase fur- 
ther (Fig. [l]^). This behavior indicates that the sharp 
transition for x = gives way to a crossover in disor- 
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FIG. 2: (Color online) Solid symbols: H — T phase boundary 
of sec deduced from specific heat data. Open symbols: field 
dependence of crossover temperatures in SCX for x = 0.0 land 
X = 0.035. The latter are defined as positions of broad specific 
heat maxima in constant-field scans. Lines are guides for the 
eye. 



dered samples, presumably to a state with short range 
order. Indeed, preliminary neutron diffraction experi- 
ments on deuterated samples, performed at the TASP 
3-axis spectrometer at Paul Scherrer Institut and at the 
D23 diffractometer at ILL, failed to observe any magnetic 
diffraction peaks associated with helimagnetic order. A 
thorough search for such reflections was performed in the 
(/i, 0, /) reciprocal-space plane for 0.45 < / < 0.55, at 
T = 50 mK and = 6 T.l^ While it can not be fully 
excluded that the peak absence is not due to disorder, 
but to a drastic change of the propagation vector along 
the chain axis, due to the large spin wave velocity this 
scenario appears unlikely. 

If the high-field phase is indeed short-range ordered, 
one can expect diffuse magnetic scattering in the form of 
anisotropically broadened peaks at Bragg positions. In 
SCX no such scattering could be found, but this result 
can not be considered conclusive. Unfortunately, experi- 
mental limitation would have made it unobservable even 
if it were present. In the disorder- free material where the 
peaks are sharp, they are extremely weak and correspond 
to only 0.04 ji^ ordered moment.^ Even a small broaden- 
ing due to short-range order in SCX would render them 
invisible above the rather high intrinsic background. 

Returning to the specific heat data in(Fig. [l] we note 
that even as the peaks broaden upon Br substitution, 
their positions are not subject to a substantial shift. The 
corresponding crossover temperatures at each field were 
identified with the maxima in C(T), obtained in empir- 
ical narrow-range Lorentzian fits. The thus determined 
crossover lines, for x = 1% and x = 3.5%, are shown in 
open symbols in Fig. |2] 

The effect of chemically introduced randomness on the 



magnetic properties is no less dramatic. The field de- 
pendence of magnetic susceptibility for the disordered 
samples are shown in Fig. [3] At all temperatures, the 
sharp peak seen for x = 0, is replaced by a broad feature 
at X = 1%, and altogether absent for x = 3.5%. This 
behavior is accompanied by an overall broadening of the 
step in the x(^) curve with increasing x. The jump in 
the temperature dependence of the magnetization is also 
smeared out by disorder (inset in Fig. [3|. Thus, magne- 
tometric measurements confirm that the BEC transition 
in the parent compound is replaced by a broad crossover 
in the disordered systems. 



III. DISCUSSION 

The effect of bond disorder on the field-induced BEC 
of magnons has been previously studied experimentally 
in a number of prototype materials, including IPA- 
Cu(Cli_^Br^)3,'^Tli_^K^CuCl3 (Ref.|20|) and most re- 
cently Ni(Cli_^Br^)2-4SC(NH2)2.^^ It was found that 
disorder may considerably affect the phase boundary, 
modify the critical behavior and give rise to a compress- 
ible disordered state just below the 3D ordering transi- 
tion. However, in all of these materials, the transition it- 
self is preserved, even for a large concentration of disorder 
sites that significantly alters the critical field and temper- 
ature. The strength and sharpness of anomalies in ther- 
modynamic properties remains practically unaffected by 
disorder. The behavior of SCX is clearly different. Here 
the condensate appears to be extraordinarily fragile with 
respect to randomness, resulting is qualitative contract 
between thermodynamic anomalies in the pure and dis- 
ordered samples. It is this contrast that we emphasize as 
the main experimental finding, since weak or broadened 
anomalies per se do not necessarily signify the absence of 
a phase transition (see, for example, Refs. 121(22)) . Below 
we propose a qualitative model that attributes this be- 
havior to a geometric frustration of magnetic interactions 
and the incommensurate nature of spin correlations. 

Consider a single geometrically frustrated spin tube 
in sec. In a magnetic field that exceeds i^c, it is a 
gapless Luttinger spin liquid^^ with a correlation length 
that diverges at T ^ 0. These increasingly long-range 
correlations will be incommensurate and helimagnetic. 
Their pitch = Qoc will be defined by the magni- 
tude of geometric frustration of interactions in the spin 
tube. A random distribution of Br sites in SCX will re- 
sult in a random distribution of the frustration ratios 
along the spin tube, and, in turn, a randomized phase 
of the spiral spin correlations. For site j, we can write 
= (j)oj + Xli<j ^^i^ where 6(j)i are defect-induced 
phase slips. Now consider the interaction of two adja- 
cent spin tubes, assuming no correlations between their 
substitution sites. Due to the random phase slips in the 
spin spirals in each spin tube, inter-tube exchange in- 
teractions will average to zero at the mean field level. 
Thus, 3D interactions that are critical for establishing 
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FIG. 3: (Color online) Field dependence of DC magnetic 
susceptibility measured in SCX with x = 0, x = 0.01 and 
X = 0.035, at T = 500 mK (a) and T = 800 mK (b). Lines 
are guides for the eye. Inset: magnetization as a function of 
temperature measured a a constant field if = 6 T > He. 
The X = data set (circles) has been shifted down by 
1 • IO^^/xb/ Cu for visualization purposes. The arrows indicate 
the positions of C/T maxima for the respective Br concentra- 
tions. 

long-range magnetic order at a finite temperature will 
be entirely wiped out by the random frustration intro- 
duced by bond disorder. This scenario is schematically 
illustrated in Fig. |4j While it is not directly supported 
by the experimental data, it appears very plausible in the 
particular case of SCX. Indeed, in our material, the inter- 
chain interactions are known to be weak.^ At the same 
time, due to the very large spin velocity (temperature 
equivalent of about 150 K), one-dimensional correlations 
are already very well developed at the experimental tem- 
peratures of a few Kelvin. In other words, upon cooling 
in a field, individual spin tubes are practically ordered, 
and it is the the inter-tube ordering that occurs at the 
transition point in the disorder-free material. 



It is important to emphasize that in the proposed sce- 
nario even weak and sparse defects will suppress 3D or- 
der. Indeed, even small phase slips induced by each im- 
purity will, over large distances, lead to an accumulation 
of phase difference between adjacent spin tubes. This is 
a specific consequence of the incommensurate and quasi- 
ID nature of the disorder- free system. 
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FIG. 4: A cartoon representation of the random frustration 
mechanism. The plot in the center shows the exchange energy 
(scalar product) of interactings spins from two adjacent one- 
dimensional helimagnets. Defects (shaded ellipses) induce 
randomly positioned phase slips, which ensures that inter- 
chain interactions average to zero. 

Incommensurability and frustration being the key in- 
gredients of our model makes it distinct from the basic 
mechanism of Bose Glasses. The latter emerges from 
the mapping the magnetic system onto a lattice gas of 
magnons with a spatially random chemical potential. 
Neither incommensurability or frustration are explicitly 
present in the resulting bosonic Hamiltonian, and are 
therefore overlooked by the approach. For our newly pro- 
posed mechanism, there are numerous theoretical ques- 
tions to be further addressed. Is there a critical strength 
of frustration disorder that is required to disrupt the 
magnon condensate? What are the critical properties of 
the order /disorder transition induced by random frustra- 
tion? How does the random frustration physics interplay 
with that of the "conventional" Bose glass? 
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